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Vaporized and multifragmenting sources produced in heavy ion collisions at inter- 
mediate energies are good candidates to investigate the phase diagram of nuclear 
matter. The properties of highly excited nuclear sources which undergo a simulta- 
neous disassembly into particles are found to sign the presence of a gas phase. For 
heavy nuclear sources produced in the Fermi energy domain, which undergo a si- 
multaneous disassembly into particles and fragments, a fossil signal (fragment size 
correlations) reveals the origin of multifragmentationispinodal instabilities which 
develop in the unstable coexistence region of the phase diagram of nuclear mat- 
ter. Studies of fluctuations give a direct signature of a first order phase transition 
through measurements of a negative microcanonical heat capacity. 



1 Introduction 

The decay of highly excited nuclear systems through a simultaneous disassem- 
bly into fragments and particles, what we call multifragmentation, is a subject 
of great interest in nuclear physics. Indeed multifragmentation should be re- 
lated to subcritical and/or critical phenomena. Thus it is fully connected to 
the nature of the phase transition which is expected of the liquid-gas type 
due to the specific form of the nucleon-nucleon interaction; as van der Waals 
forces, the nucleon-nucleon interaction is characterized by attraction at long 
and intermediate range and repulsion at short range. 

Although multifragmentation has been observed for many years, its ex- 
perimental knowledge was strongly improved only recently with the advent of 
powerful devices built in the last decade. Selecting the "simplest" experimen- 
tal situations, well defined systems or subsystems which undergo vaporization 
(simultaneous disassembly into particles) or multifragmentation can be thus 
identified and studied. 

It is a difficult task to deduce information on the phase diagram and the 
related equation of state of nuclear matter from nucleus- nucleus collisions at 
intermediate energies. But it is also a very exciting novel physics in relation 
with thermodynamics of finite systema^i^connection to other fields) without 
external constraints (pressure,volume) Era. 
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2 From multifragmentation to vaporization : identification of 
the gas phase 
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Figure 1. Fragment mul- 
tiplicity (normalised to the 
number of incident nucle- 
ons) as a function of theiex- 
citation energy, (from H). 
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Let us first locate in which excitation energy domain multifragmentation 
takes place. Fig ^ indicates the evolution of the reduced (normalised to the 
size of the multifragmenting system) fragment multiplicity as a function of 
the excitation energy per nucleon of the system. A universal behavior char- 
acterized by a bell shape curve is observed. The onset of multifragmentation 
is observed for excitation energies around 3 MeV per nucleon, the maximum 
for fragment production is found around 9 MeV per nucleon, i.e. close to 
the binding energy of nuclei. At higher excitation energy, the opening of the 
vaporization channel reduces fragment production. 
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Figure 2. Composi- 
tion of vaporized quasi- 
projectiles, formed in 95 
AMeV '^^Ar+^^Ni colli- 
sions, as a function of their 
excitation energy per nu- 
cleon. Symbols are for data 
while the lines (dashed for 
He isotopes) are the re- 
sults of the model. The 
temperature values used in 
the m«dcl are also given, 
(from □). 



The gas phase has been identified by studying the deexcitation properties 
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of vaporized quasi-projectiles with A around 36 Chemical composition (first 
and second moments) and average kinetic energies of the different particles 
are well described by a gas of fermions and bosons in thermal and chemical 
equilibrium. Inclusion of a van der Waals-like behavior (final state excluded 
volume) was found decisive to obtain the observed agreement (see for example 
figure 1^). In the model, the experimental range in excitation energy per 
nucleon of the source was covered by varying the temperature from 10 to 25 
MeV and the free volume was fixed at 31^0: which corresponds to an average 
inter-distance between particles of about 2 fm, close to the range of the nuclear 
force (freeze-out configuration). 

3 Thermometry and calorimetry : caloric curves and 
first-order phase transition ? 

The plateau observed in the shape of the caloric curve (determined from 
calorimetry and nuclear thermometry) was proposed by the ALADIN-collab- 
oration a few years ago as a signature of a first-order phase transition Since 
this observation works from different cnllahorations, covering a large range in 
mass of nuclei, have been published □'Lracl. Many caloric curves have been 
obtained which can roughly be classified in two groups depending on the nu- 
clear thermometer chosen (isotopic double ratios using Li/'^ Li) / (^He/^He) 
or [d/i) / He/'^He)). Moreover the presence of a plateau was not confirmed 
in these studies, even by the ALADIN collaboration when looking, at proper- 
ties of target-like spectators in Au-f Au collisions at lOOOAMeVEd. All these 
studies clearly indicate that no decisive signal can be extracted. We do not 
have an absolute nuclear thermometer and above all, experimentally one does 
not explore the caloric curve at constant pressure nor at constant volume. In 
fact measured caloric curves are sampling a monodimensional curve oiL±hfi 
microcanonical equation of state surface (T versus energy and volume) BEj; 
for each energy of the system a different average volume at freeze-out (no 
more nuclear interaction) is obtained, depending on the observed partition. 

4 Statistical and dynamical descriptions of multifragmentation 

Many theories^ have been developed to explain multifragmentation (see for 
example ref. t3 for a general review^ models). Among the models some are 
related to statistical approaches Ej'Eil, valid at and after freeze-out, whereas 
others try to describe the dynamical evolution of systems, from thfLJaeginning 
of the collision between two nuclei to the fragment formation li3'E£l. I shall 
very briefly focus here on two models which will be compared later on to ex- 
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perimeptal data. Firstly a statistical description of multifragmentation (SMM 
model llil), in which an equilibration of a system at low density is assumed. 
Then the statistical weight of a given break-up channel /, i.e. the number 
of microscopic states leading to this partition, is determined by its entropy 
AT f—exp S f within the microcanonical framework. In such an approach the 
initial parameters as the mass and charge of the multifragmcnting system, 
its excitation energy, its volume (or density) and the eventual added radial 
expansion have to be backtraced to experimental data. Secondly, dynamical 
stochastic mean-field simulations which are obtained by restoring fluctuations 
in deterministic one-body kinetic simulations. In particular in such simula- 
tions, relative to the standard nuclear Boltzmann treatment, an approximate 
tool is provided by introducing a noise by means of.aJprownian force in the 
mean field (Brownian One-Body (BOB) dynamics tHO). The magnitude of 
the force is adjusted to produce_Lhe same growth rate of fluctuations as the 
full Boltzmann-Langevin theory 13. Such simulations completely describe the 
time evolution of the collision and thus help in learning about nuclear matter 
and its phase diagram whereas statistical models start from the phase diagram 
and have more to do with the thermodynamical description of finite nuclear 
systems. 

Both descriptions have been successful in reproducing peerage static and 
kinematical properties of fragments (see for example ref. c3). They will be 
compared in what follows to more constrained observables which are expected 
to bring decisive information on the origin and properties of multifragmenta- 
tion. 

5 Correlations in events: spinodal instabilities and 
equilibration 

5. 1 Fragment size correlations: a fossil signature 

Dynamical simulations predict that during a central collision between heavy 
nuclei in the Fermi energy domain (30-40 MeV per nuclcon incident energies) 
a wide zone of the phase diagram is explored (gentle compression-expansion 
cycle) and the fused system enters the liquid-gas coexistence region (at low 
density) and even more precisely the unstable spinodal region (domain of 
negative incompressibility). Thus a possible origin of multifragmentation may 
be found through the growth of density fluctuations in this unstable region. 
Within this theoretical scenario a breakup into nearly equal-sized "primitive" 
fragments should be favored in relation with the^ wave-lengths of the most 
unstable modes present in the spinodal region til. However this picture is 
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expected to be blurred by several effects: the beating of differents modes, the 
presence of large wave-length instabilities and eventual coalescence of nascent 
fragments. Then how to search for a possible "fossil" signature of spinodal 
decomposition? A few years_ago a new method called higher order charge 
correlations was proposed in l3. All fragments in one event (average fragment 
charge < Z > and the standard deviation per event AZ) are used to build 
the charge correlation for each fragment multiplicity. 



Experiment BoB 




AZ 




Figure 3. Fragment 
charge correlations from 
fused events produced in 
central collisions between 
Xe and Sn at 32 MeV per 
nucleon incident energy: 
comparison between exper- 
iment (left) and BOB cal- 
culations (right) for frag- 
ment multipliciti«s equal to 
4 and 6. (from £3). 



Figure |^ shows results from E2I for such correlation functions in experi- 
mental fusion events and BOB simulated events (Xe-|-Sn system at 32 MeV 
per nucleon). For all fragment multiplicities the charge correlation has a peak 
in the bin AZ — 0-1, indicating an enhancement of partitions with equal- 
sized fragments. This weak but non ambiguous enhancement (0.1% of events 
if we restrict to the bin 0-1 and about 1% if we enlarge to bin 1-2 to take 
into account secondary decays of fragments) is interpreted as a signature of 
spinodal instabilities as the origin of multifragmentation in the Fermi energy 
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domain. Moreover the occurrence of spinodal decomposition signs the pres- 
ence of a hquid-gas coexistence region and consequently, although indirectly, 
a first order phase transition. 

5.2 Fragment-particle correlations: equilibrium at freeze- out 

Fragment-particle velocity correlations in events have been proposed to ex- 
perimentally measure excitation energy of hot primary fragments produced 
in multifragmentation c3. By means of this technique multiplicities and rela- 
tive kinetic energy distributions between fragments and light charged particles 
that they evaporate are determined to reconstruct the excitation energies of 
fragments. For the Xe-|-Sn reaction, the INDRA collaboration has measured 
the evolution of the average excitation energy per nucleon of primary frag- 
ments produced in multifragmentation o£_£used systems at different incident 
energies (from 32 to 50 MeV per nucleon) LJ. Within the error bars a constant 
value around 3.0-3.5 MeV per nucleon was measured in good agreement with 
the approach at equilibrium (SMM). This suggests that equilibrium is reached 
at freeze-out. Note that dynamical simulations (BOB) performed at 32|_MeV 
per nucleon also predict the same excitation of fragments at freeze-out cB. 

How to reconcile the dynamical (spinodal instabilities) and sta- 
tistical (equilibrium at freeze-out) aspects which have been ex- 
tracted from correlations ? The following scenario can be proposed: 
spinodal instabilities cause multifragmentation but when the system 
reaches the freeze-out stage, it has explored enough of the phase 
space in order to be describable through an equilibrium approach. 

6 Kinetic energy fluctuations and negative microcanonical 
heat capacity 



Within the microcanonical equilibrium framework, it was recently shown 




that for a given total energy of a system, the average partial energy stored in a 
part of the system is a good microcanonical thermometer, while the associated 
fluctuations can be used to construct the heat capacity (see □) . In presence of a 
phase transition large fluctuations are expected to appear as a consequence of 
the divergence and of the possible negative branch of the microcanonical heat 
capacity. From experiments the most simple decomposition of the total energy 
E* is in a kinetic part Ei and a potential part E2 (Coulomb energy -f total 
mass excess). However these quantities have to be determined at freeze-out 
and consequently it is necessary to trace back this configuration on an event 
by event basis. The true configuration needs the knowledge of all the charged 
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Figure 4. Measurements of microcanoni- 
cal heat capacity per nucleon (symbols) as a 
function of the excitation energy per nucleon 
for quasi-projectiles produced in Au+Au coUi- 
sions. The two panels refer to different freeze- 
out hypotheses. The grey contouiHpdicates 
the confidence region for Ct- (from 
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Figure 5. Same as figure ^for fused nuclei 
produced in central collisions between Xe 
and Sn at 32 rMeV per nucleon incident 
energy, (from E3). 



particles evaporated from primary hot fragments and of the undetected neu- 
trons; consequently some reasonable hypotheses have to be done. Note also 
that fragment-particle correlations discussed just before can help to obtain a 
better knowledge of freeze-out configurations (see c2i) . Then the experimental 
correlation between the kinetic energy per nucleon Ei/Aq and the total ex- 
citation energy per nucleon E* /Aq of the considered system can be obtained 
as well as the reduced variance of the kinetic energy af/ < Ef >. Finally 
the microcanonical temperature of the system can be obtained by inverting 
the kinetic equation of state and the total microcanonical heat capacity Ct is 
extracted from the following equations: 

and C- 



5T " c 



Figures ^ and || show results obtained by M. D'Agostino et al. and the 
INDRA collaboration for hot nuclei with mass number around 200 formed in 
different experimental conditions. In figure ^the micocanonical heat capacity 
is calculated over a large excitation energy range for quasi-projectiles (formed 
in Au-|-Au collisions at 35 MeV per nucleon incident energy) assuming two 
different hypotheses to trace back freeze-out configurations; figure || refers to 
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fusion events produced in central Xe+Sn collisions at 32 MeV per nucleon; in 
this latter case, a narrower excitation energy distribution (bell shape curve on 
the figure) is observed. A distinct negative branch is observed, revealing a first 
order phase transition. The distances between the poles are associated with 
the latent heat. Note that the same location of the pole at high excitation 
energy is found when similar hypotheses are made for freeze-out reconstruction 
(left part of figure ^ and figure ||). 

7 Signatures of critical behavior 

For finite systems, related to the correlation length, a critical region instead 
of a critical point is expected. The following signatures of critical behavior 
were reported: power laws have been observed within selected conditions, 
critical exponents—have been measured in agreement with with those of a 
liquid-gas model tBc3 assuming that fragment multiplicity or thermal excita- 
tion energy is the control ( "order" ) parameter pad a nuclear scaling function 
has been evidenced by the EOS collaboration c3. However, as compared to 
infinite systems, potential divergences are smoothened over finite regions of 
the chosen control parameter and the choice of the fit regions where the data 
are assumed to reflect the critical behavior is crucial. Clearly one needs a 
precise and objective procedure in order to determine order parameters and 
critical signals. Such a methodology was recently proposed for second order 
phase transitions EJ and we can expect in the future to dispose of a similar 
methodology for first order phase transitions. 

8 Conclusions 

A set of coherent results showing the existence of a first order phase transition 
in nuclear matter has been obtained and the two signals observed related to 
correlations and fiuctuations constitute a strong starting point for systematic 
investigations. Clearly caloric curves do not and can not give a decisive sig- 
nal of a first order phase transition. Selected data have properties compatible 
with the equilibrium hypothesis at freeze-out and this framework is up to now 
chosen for progressing in the experiment-theory interaction. Experimentally 
an effort has to be made to better define configurations at freeze-out, which 
are key points to bring more quantitative information (latent heat. . . ). On 
the theoretical side, concerning thermodynamics of finite systems, some im- 
provements are needed in lattice gas models for example to take into account 
the specificities of nuclei (quantal aspects and Coulomb interaction). Con- 
cerning the signatures of critical behavior and the definition of the critical 
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region, experimentalists need a methodology dedicated to first order phase 
transition for finite systems. 

/ am highly indebted to my colleagues of the INDRA collaboration and to 
R. Botet, X. Campi, Ph. Chomaz, M. Colonna, M. D'Agostino and F. Gul- 
minelli for valuable discussions. 
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